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Abstract 



' The reaction pn d(f> is studied within a covariant boson exchange model. 

' The behavior of polarization observables being accessible in forthcoming ex- 

. periments near threshold is predicted. 



. PACS. 13.75.-n Hadron-induced low- and intermediate-energy reactions and 

scattering (energy < 10 GeV) - 14.20.-c Baryons (including antiparticles) - 
21.45.+V Few-body systems 
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I. INTRODUCTION 



> 
(N 

^ I Data on elementary reactions with neutrons are scare since either they must be extracted, 

Q I with some efforts and even mostly with some model dependent assumptions, from reactions 



on nuclei, or a tagged neutron beam (cf. |]T|) is used. The spectator technique represents 
one example how one can use a deuteron target to isolate quasi-free reactions at the neutron. 
It is based on the idea to measure the spectator proton (psp) at fixed (or slightly varying) 
beam energy in the meson (M) production reactions pd —>■ dMpsp thus exploiting the internal 
^ I momentum spread of the neutron inside the deuteron {d). In such a way one gets access 
C ■ to quasi-free reactions pn — *• dM if, in experiments with the deuteron target at rest, the 

spectator proton has momenta in the order of 50 ■ ■ ■ 150 MeV/c. 
^ , An experimental investigation of the threshold-near (pseudo) scalar and vector meson 

^ ! production at the neutron becomes therefore feasible. Indeed, at COSY the ANKE spec- 
trometer set-up can be used, in particular, for studying the cq, lu and production with the 
internal beam at a "neutron target" 0. This offers the possibility to enlarge the data base 
on hadronic reactions and to address special issues. For instance, there is already a large 
body of data which can be used for a systematic study of the OZI rule violation via u and (j) 
production in vrA^ and pp reactions (cf. 0] for a reanalysis) and in pp annihilation (cf. 
for theoretical analyzes) as well. OZI rule violations are of interest with respect to possible 
hints to a significant ss admixture in the proton, as supported by the pion-nucleon S term 
and interpretations of the lepton deep- inelastic scattering data . Besides the impact 



on hadron phenomenology the origin of the OZI rule has also a link to QCD |]TQ,11 . 



*0n leave of absence from Bogoliubov Laboratory of Theoretical Physics, JINR, 141980, Dubna, 
Moscow reg., Russia 
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In the reactions pp ppM or pn pnM the final state interaction among the outgoing 
nucleons plays a role. Therefore, the meson production process is a convolution of the "pure 
production process" and the final state interaction. In the case of the reaction pn dM 
one has one well defined final state of the nucleons and may better constrain the elementary 
production amplitude. 

Finally, we mention that in the reactions pp —>■ ppM, pn —>■ pnM and pn —>■ dM the 
conservation laws and symmetry principles determine a different dynamics near threshold, 
which needs to be investigated to allow a firm understanding of the reactions and the sys- 
tematics of the OZI rule violation. 

Given this motivation, in [0,0] the reaction pn dV with V = uj,(f) has been studied 
in some detail, {pn dS with S = aQ,-!],-)]' is considered in [|14|-) In [0] the cross sections 



and angular distributions are elaborated as a function of the excess energy within a two-step 
model. The same observables are evaluated in |jl3|] within the framework of a boson exchange 
model with emphasis on the ratio of cross sections <Jpn-^d4>/o'pn-^du} being of direct relevance 
for the OZI rule violation. 

Since at COSY the proton beam is polarized and also polarized targets are envisaged, 
we extend the previous studies [p!2| , p!3[ to make a prediction of polarization observables in 
the reaction pn —>■ dcj). The set-up of the ANKE experiment at COSY can directly identify 
the (f) via its K^K~ decay channel. We present the asymmetry, tensor analyzing power, 
and proton-phi spin-spin correlations. In doing so we use our previous one-boson exchange 
model |]15| with parameters adjusted to available threshold-near data on production in pp 
and irp reactions and combine this with our previous studies |1T5|-|T7|] employing a solution 
of the deuteron wave function within the Bethe-Salpeter (BS) formalism. In this way we 
derive a completely covariant approach with correct treatment of the off-shell effects in the 
vrA^ — > (j)N subprocess amplitudes. 

Extensions to pn — > du and even to pd — > dVpsp are much more involved and, therefore, 
relegated to separate work. 

Our paper is organized as follows. In section 2 we present the theoretical framework and 
elaborate the basic equations. The numerical results and their discussion are presented in 
section 3. The summary can be found in section 4. 



II. THE MODEL 

The invariant differential cross section of the reaction pn d(f) reads 

where s is the square of the total energy of the colliding particles p and n in the center of 
mass, t is the square of the transfered 4-momentum, m is the nucleonic mass, Si, S2, M.^ and 
M-d denote the spin projections on a given quantization axis, and T stands for the invariant 
amplitude. The general form of T may be written in the form (cf. Fig. T^st^'^{s^t) = 
{ d, Aid \Gii 2 ), where is the polarization 4-vector of the meson. The scattering 

operator G represents a 4-vector in Minkowski space, a vector in the vector space of mesons 
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and a 16® 16 component object in the spinor space of nucleons. The deuteron is described as 
a 16 component BS amphtude $(1,2) which is defined as a matrix element of a time ordered 
product of two nucleon fields i){x) by $"^(1, 2) = ( d\T (z^"(l)V^^(2)) |0 ) and satisfies the BS 

equation. By defining another scattering operator via O = the invariant amplitude 

reads 

T^st^^{s,t) = -^ J ^$^^(1', 20 0^^,(12; 1%^,)<(1)<(2), (2) 

where $^^(1',2') is the conjugate BS amplitude in the momentum space, p is the relative 
4-momentum of the nucleons in the deuteron, and u{l) and u{2) denote the Dirac spinors 
for the incident nucleons. Summation over spin indices a, f5,b,c = 1 ■ ■ • 4 occurring pairwise 
is supposed. The operator 0^*^(12; 1'2', A4(f,) is a scattering operator describing the meson 
production in the final state. This operator acts in the spinor space of protons and neutrons 
separately; the upper (Latin) and lower (Greek) spinor indices refer to protons and neutrons, 
respectively. The first indices, b and a, form an outer product of two columns, whereas the 
second ones, c and j3, form an outer product of two rows. To specify explicitly the spinor 
structure we decompose the operator O in each of its indices over the corresponding complete 
set of Dirac spinors, i.e, 

d^^^(12; 1'2',M^) = ^ E ^;:'':pp'(12; l'2')«^(l')n^(l)n?(2)<,(2'), (3) 



r,r' ,p,p'=l 



where the coefficients A^,*^ ,(12; 1'2') may be found by using the completeness and orthog- 



rr' ,pp' 



onality of the Dirac spinors, Ur{p)ur'{p) = 2erm6rr'^ yielding 

A^,'^^^,{12; 1'2') = Srer'SpSp' 4(1') ^'(2') d%{12- 1'2' , M^) <(1) «J(2), (4) 
where = +1 for r = 1, 2 and = — 1 for r = 3, 4. Substituting (|^) into (g) one obtains 

^irr^(^>^) = 7£p/ E <t,.^(i2;i'2')<,(2')«::(i') 



(2m)2y (27r) 
2 " d^p 



(2m 



,{12-I'2')ul{l') 



= E / (^<:..'(12; l'20^;.K2')^A..(l',2')n.(l'), (5) 

where 7c is the charge conjugation matrix, Vr{2') = (Mr(2'))"^7c, and the new BS amplitude 
^_M^{1',2') = 7c$x^(l', 2') now is a 4 (g) 4 matrix and represents the solution of the BS 
equation written also in matrix form. 

For further evaluations of the amplitude (|^) one needs to specify an explicit form of the 
operator O. Following |]l5l, and taking into account only the meson-exchange contribution 
from one-boson exchange diagrams with internal meson conversion, the operator O may be 
represented by four truncated diagrams as depicted in Fig. |^. (The contributions of the 
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nucleonic currents are found to be negligibly small [T^JT^ and, therefore, can be omitted 



for many purposes.) Consequently, having computed the operator O by using the effective 
interaction Lagrangians for the ttNN, pNN, (ppir vertices and monopole form factors for 
dressing the vertices as in [Q, it is straightforward to obtain the coefficients 74^/^^,(12; 1'2') 

in (H). If all particles were on mass shell, A^,'^pp, (12; 1'2') exactly coincides with the amplitude 
of the elementary process 1 + 2 —>■ 1' + 2' + (p. However, in general case this amplitude 
corresponds to a virtual process of vector meson production with two off-shell nucleons in 
the final state. 



Since our numerical solution |T^ of the BS equation has been obtained in the deuteron's 
center of mass, all further calculations will be performed in this system. First, as depicted 
in Fig. |I|, we introduce the relevant kinematical variables as follows: pi^2 are the 4-momenta 
of incoming nucleons, p'l ^ stand for the 4-momenta of the internal (off-shell) nucleons in the 
deuteron with p = (p'^^ — p2)/2; denotes the polarization 4- vector of the deuteron. In 



this notation the BS amplitudes in the deuteron's rest system are of the form |17 



*iJ,^W,P'2) =Ar(fci + m)l±^U.(^2-m)05(po,|p|), (6) 

*:M,^(p'i>P2) = "^(^1 + {^^i + 2^*-^^ ~ ^2)(p6f)^ (^2 - m)0B(PO, IpI), 

where " means contraction with Dirac matrices, and k\^2 are on-shell 4-vectors related to the 
off-shell vectors p'^ 2 follows 



fci = (^p,p), k2 = {Ep, -p), p[ = {p[q,p), P2 = ip'20, -P), Ep = ^p'^ + m?, (7) 

and 0s,d(po, |p|) are the partial scalar amplitudes, related to the corresponding partial ver- 
tices as 

, / I IN Gs,d{po, IpI) , . 

(f^sAvo, IpI) = 72 • (8) 



(^^rf - Ep 



Md is the deuteron mass, and the normalization factor is A/" = ^'\/^2E[E + m)| . To be 
explicit let us recall the components of the polarization 4-vector of a vector particle with 
4-momentum p = {E,p), polarization index Ai = ±1, and mass M as 

= + ^M(?+M)) ' 

where ^j^^ is the polarization 3-vector for the particle at rest with = ^(l,i,0), = 
■^(1,— z,0), £q = (0,0,1) and the above Dirac spinors, normalized as u{p)u{p) = 2m and 
v{p)v{p) = —2m, read 



u{p, s) = yM + e [op , v{p,s) = Vrr^ , (10) 
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where Xs = ~'i<^yXs^ and Xs denotes the usual two-dimensional Pauli spinor. In general, the 
BS amplitude consists on eight partial components. In we take into account only the 
most important ones, namely the S and D partial amplitudes. The other six amplitudes may 
become important at high transferred momenta [l^Jl^] ? hence for the present near-threshold 
process they may be safely disregarded. Substituting (|6|-p!0|) into @ one obtains after some 
algebra 



in 



r,r' 



\Md\ + 1 



(11) 



2 :Aj2,rr'{PlP2,PlP2)(>r+r',Ma 



Po 



G 



D 



A^S2,rr'iPlP2,P'lP2)xt' ('^^)Xr (^^Xl^^ 



Po 



where here and below the sums over r, r' run from 1 to 2. By closing the integration contour 
in the upper hemisphere and picking up the residuum at pq = Md/2 — Ep and introducing 
the notion of the deuteron 5* and D wave functions as 



usip) 



Gs{po, \p\) 



ud{p) 



Gd{Po, \p\) 



47rvOTrf(2Ep - Ma) ' '"^^ 47rvOTd(2Ep - M^j 
with 2 J d\p\ + ^ IT, the final expression for the amplitude may be written as 

d^p 



(12) 



SlS2 



S,t) 



\Md\ + l 



us{p) - 



Epi2'n 
ud{p) 



V2 



^^siS2,rr'iPlyP2'j Pi ~P 

ud{p) 



r+r',Md 



3^ (eMd^^) xt'{^n)xA , (13) 



where n is a unit vector parallel to p. 



III. DISCUSSION OF RESULTS 

In our evaluation of the above equations we use, for describing the deuteron wave function 
, the numerical solution of the BS equation in ladder approximation obtained with a 
realistic one-boson exchange interaction which includes 7T,a,uj,6, p and rj exchanges. The 
effective parameters used in the ladder approximation have been fixed in such a way to 
obtain a good description of the NN elastic scattering data and the static properties of 
the deuteron ||T^. Independent of this set of parameters related to deuteron properties, 
the effective coupling constants and the attributed cut-off factors are taken from the recent 
analysis [0. In all subsequent numerical calculations, the set B from is used. For 
this parameter set the meson exchange term is by far dominating. In Fig. |^ the total cross 
section near the threshold is depicted. The shape of our cross section is rather similar to 



the one computed in [13|. However, there is a difference in the absolute values by roughly 



a factor of 2/3. This difference is in the same order of magnitude as the difference of the 
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previously often used cross section a = 0.26 /ib deduced from and the published value 
a = 0.19 fih ||20] for the reaction pp ppcf) at an excess energy of 83 MeV. 



Refs. |13|,|T5| show that several sets of parameters equally well describe the pp — 
data. These sets differ not only by absolute values of parameters but also by the relative 
contributions of meson current and the nucleon current terms. Since in case of the pn 
d(j) processes the isospin transition corresponds to A/ = the meson-exchange diagrams 
are enhanced by a factor of three in comparison with the nucleon current terms. This 
means that (i) the contributions of the nucleon current are suppressed by about one order 
of magnitude in comparison with the meson current, and (ii) the behavior of the cross 
section and angular distribution in the process pn — > dcp is expected to follow essentially the 
behavior of the dominating meson-exchange current contribution in the elementary processes 
pn — >■ pn(j). The occurrence of the deuteron wave function will only modify this behavior. 
This is clearly seen in Fig. ^, where the behavior of the angular distribution is very similar to 
the distribution found in |T3] for the reaction pn pncj). At the threshold the distribution 
is fairly flat, while with increasing excess energy some forward-backward peaking becomes 
visible. This feature depends upon the parameter set used. For instance, if we were using the 



set C of 1 15], then a slight suppression of the cross section would be expected in the forward- 



backward directions as predicted for the elementary cross section ||T^. This sensitivity of the 
predicted angular distribution has been pointed out in [13,21| as a tool to constrain further 
the parameters, if precise data become available. 

Let us now discuss the polarization observables. Near the threshold, in the final state one 
has: J/ = 0, L/ = 0, 2, JJ = + J]^ = 0~, 1~, 2~, where /, L and .P are the total isospin, 
total radial angular momentum, angular momentum and parity, respectively. Thus from 
symmetry constraints, in the initial state the allowed configurations are Jj = 0, Lj = 1, 3, 5 • • • 
and total spin Si = 0. The conservation law implies L, = 1, so that J[ = 1~ . From these 
transitions in \T^st''^'^{s,t)\'^ one may form different combinations of spin observables, which 
near the threshold behave quite differently. If one considers, for example, the cross section 
averaged over all final projections A4 at different initial spin projections, then the beam- 
target asymmetry, defined here as 

^ _ CT(Sp + Sn = l) + (rjSp + Sn = -1) - Cr{Sp + Sn = 0) 
0"(Sp + S„ = 1) + a{Sp + Sn = -1) + Cr(Sp + Sn = 0)' 

is predicted to be —1 near threshold and to increase with increasing energy. This is in 
contrast with the asymmetry defined for the elementary process pp — > pp0 which is -|-1, 
as predicted in |]T5[, and decreases with increasing energy. The energy dependence of the 
asymmetry (|l^ is depicted in Fig. |^. 

Another interesting polarization observable is the tensor analyzing power, which may be 
defined either for the final deuteron or for the final meson. For instance, the deuteron tensor 
analyzing power is defined by the cross section averaged over all spins but the deuteron, 

1 a{Md = I) + (T{Md = -1) - 2a{Md = 0) . . 

'° ^ a{Md = l) + cr{Md = -l) + a{Md = 0)- ^ ^ 

It is seen from ( |T5| ) that in line with the above discussed selection rules, cr(A^rf = 0) -C 
<j{M.d = ±1) and the deuteron tensor analyzing power is predicted to be almost constant, 
T20 ~ l/"\/2, in a large region of the energy excess. This prediction is illustrated in Fig. |^. 
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In Fig. 1^ another polarization observable is depicted defined by 



W, 



zz 




1) 



Sp=+l/2 



(16) 



which characterizes the proton-meson spin-spin correlation. This quantity is predicted to 
vanish near the threshold for the meson exchange diagrams. A substantial deviation of Wzz 
from zero would point to the presence of non-negligible nucleon current contributions in the 
production process. In Fig. ^ it is seen that there is some weak dependence of W^z upon 
the energy excess. Nevertheless Wzz remains very small indicating that the polarizations 
of the incident proton and outgoing meson are almost uncorrelated. A different situation 
will occur in the reaction pd —>■ d(j)psp. In this case, if one measures also the polarization 
of the spectator proton, the quantity Wzz is expected to strongly depend on whether the 
polarizations of protons are the same or have opposite directions. 

Since we have T^st''^''{s,t) at disposal, via (|l]) any other polarization observable is ac- 
cessible within the our formalism and corresponding numerical code. 



In summary we present here a evaluation of polarization observables for the process 
pn — i> d(j) which are accessible in forthcoming experiments. Off-shell effects in the sub- 
process vrA^ (f)N can correctly be dealt with if one restricts the treatment on the by 
far dominating meson-exchange current. It is the beam-target asymmetry which differs 
drastically from the one in the reaction pp — > pp(j). The tensor analyzing power is fairly 
insensitive to variations of the excess energy, and the proton-phi spin-spin correlation is 
very small. 

The treatment can be extended to incorporate the u meson production; work along this 
line is in progress. 
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FIGURES 




FIG. 1. The diagram for the process p{l) + n(2) = d + (p. G is the scattering operator, and ^* 
denotes the polarization vector of the outgoing vector meson. 




FIG. 2. Graphical representation of the operator O defined in (^) in the one-boson exchange 
approximation. Nucleons in the initial and final states are represented by thin lines which are 
truncated by vertical dotted lines as to obtain an operator. The exchanged vr and p mesons 
are depicted by vertical dashed and thick lines, respectively. The four different contributions 
correspond to different combinations of and p^,p^ exchanges. 
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FIG. 3. Total cross section for the reaction pn d(f) as a function of the excess energy 
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FIG. 4. Angular distribution in the center of mass system for various values of the excess energy. 
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FIG. 5. The beam-target asymmetry ^ as a function of the the excess energy. 
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FIG. 6. Deuteron tensor analyzing power T20 as a function of the excess energy. 
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FIG. 7. Proton-phi spin-spin correlation Wzz as a function of the excess energy. 
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